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Starting from the Boltzmann equation we calculate the 
frequency and the damping of the monopole and quadrupole 
oscillations of a classical gas confined in an harmonic poten- 
tial. The collisional term is treated in the relaxation time ap- 
proximation and a gaussian ansatz is used for its evaluation. 
Our approach provides an explicit description of the transition 
between the hydrodynamic and collisionless regimes in both 
spherical and deformed traps. The predictions are compared 
with the results of a numerical simulation. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40. Db 



I. INTRODUCTION 

After the experimental realization of Bose-Einstein 
condensation in trapped atomic gases [Q], the investiga- 
tion of collective excitations in these systems has become 
a very popular subject of research (see for a recent 
theoretical review). At very low temperature, when the 
the whole system is Bose-Einstein condensed, the motion 
is described by the hydrodynamic equations of super- 
fluids. These equations, which can be directly derived 
starting from the mean field Gross-Pitaevskii equation 
for the order parameter, give predictions || in very good 
agreement with experiments . At higher temperatures 
the mean field effects become less important, while col- 
lisional terms cannot be longer ignored. If the temper- 
ature is notably larger than the critical temperature for 
Bose-Einstein condensation the dynamical behaviour of 
a dilute gas is well described by the Boltzmann equa- 
tion. In this case two different regimes may occur: a 
collisional (hydrodynamic) regime characterized by con- 
ditions of local statistical equilibrium and a collisionless 
regime where the motion is described by the single par- 
ticle hamiltonian. Differently from the case of uniform 
gases, also in the collisionless regime the system exhibits 
well defined oscillations which are driven by the exter- 
nal confinement. The equations for the collisional regime 
were investigated in (^||, while a phenomenological in- 
terpolation between the two regimes was proposed in . 

The purpose of this paper is to provide an analysis of 
the lowest oscillation modes (eigenfrequency and damp- 
ing) using the classical Boltzmann equation. In partic- 
ular we study the transition between the hydrodynamic 
and collisionless regimes. Our approach is based on the 
relaxation time approximation for the collisional integral. 



The corresponding relaxation time is explicitly calculated 
by means of a gaussian ansatz for the distibution function 
which allows for an analytic evaluation of the collisional 
contribution to the equations of motion. 

For harmonic trapping the gaussian ansatz exactly re- 
produces the solution of the classical Boltzmann equa- 
tion in both the collisional and collisionless regimes and 
is consequently expected to be a good approximation also 
in the intermediate regime. This is explicitly checked by 
comparing the analytic predictions of the gaussian ansatz 
with the exact results of a numerical simulation based on 
molecular dynamics. 

In our paper atoms behave like hardspheres, <7o being 
their total cross section which will be assumed to be en- 
ergy independent. This is well satisfied in classical ultra- 
cold gases where collisions are completely characterized 
by the s-wave scattering length and the cross section is 
thus isotropic and in most cases energy independent. 



II. METHOD OF AVERAGES 

The starting point of our analysis is the Boltz- 
mann equation for the phase space distribution function 
Z(r,v,i) §: 



df 
dt 



V v / = IcoIl(/) 



(1) 



where 



/odl(/(r,Vi)) = g : J d 2 nd 3 v 2 |v 2 -vi| 

[/(vl)/K) - /(vi)/(v 2 )] 

is the usual classical collisional integral. 

It accounts for elastic collisions between particles 1 and 
2, with initial velocities vi et V2, and final velocities 
and v 2 . The solid angle £1 gives the direction of the fi- 
nal relative velocity. The expression for the collisional 
term can be easily extended to include effects of both 
Bose and Fermi statistics [@]. Actually most of the re- 
sults discussed in this paper hold also in the presence 
of quantum degeneracy, provided the system is not Bose- 
Einstein condensed and one can ignore mean field effects. 
The quantitative estimates of collisional effects will be 
however based on classical statistics. 
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The force F tra p = — VC/trap(£, y, z) is produced by the 
confining potential which in the following will be chosen 
to be of harmonic form: 

t/trapK y, z) = ^mu}'i(x 2 + y 2 ) + ]^nuj 2 z z 2 . (2) 

We introduce the anisotropy parameter A = uj z /loj_. For 
A = 1, one deals with an isotropic harmonic trap, for 
A 2 « 1 one has a cigar shaped trap and for A 2 > 1 a 
disk shaped trap. Starting from (Q), one can derive useful 
equations for the average of a general dynamical quantity 
X(r,v): 

^ - <v • V rX ) - • V vX ) = (x/coii) , (3) 

at m 

where the average is taken in both position and velocity 
space: 

(X) = ^ J d 3 r d 3 v /(r,v,i) X (r,v) , (4) 

The quantity (xlcott) can be written in the useful form 

(X^coii) = f d 3 r d 3 Vl [xi +X2- Xv - X2'} Icott(f) , (5)' 

as a consequence of the invariance properties of the cross 
section (we have put Xi — x( r i v i))- The collisional 
contribution (^J) is equal to zero if x corresponds to a 
dynamic quantity conserved during the elastic collision. 
This happens if x can be written in the form |s|]lC|| 



X = a(r) + b(r) • v + c(r)v 2 



(6) 



force and for the energy. The collisional term does not 
contribute to the above equations because all the dy- 
namic quantities satisfy the criterium (^|). So there is no 
damping for the "breathing" mode of a classical dilute 
gas confined in an harmonic isotropic trap. The same 
is true if one includes quantum degeneracy effects in the 
collisional term of the Boltzmann equation. 

By looking for solutions of eqs.((7-^) evolving in time 
as e luJt one immediately finds the result ui = 2ojq, hold- 
ing for all collisional regimes from the collisionless to the 
hydrodynamical one. The occurrence of this monopolc 
undamped solution was first pointed out by Boltzmann 
(see, for example, the discussion in JTT[] ) . 

It is worth noticing that the frequency of the classical 
monopole oscillation in isotropic harmonic traps differs 
from the one of a Bose-Einstein condensed gas at T = 0. 
In the latter case the monopole oscillation is still un- 
damped, but the frequency, for large N, is u> — V5ojq 
||. The difference is the consequence of the mean field 
interaction which is absent in the Boltzmann equation 
(|l|). Furthermore, at finite temperature the monopole 
oscillation is expected to exhibit damping because of the 
coupling between the condensate and the thermal com- 
ponent of the gas. 



III. MONOPOLE OSCILLATION IN HARMONIC 
ISOTROPIC TRAPS 

Let us consider an harmonic isotropic trapping poten- 
tial (uj x = Ljy — lu z = luq). As a first application of 
(||), one can immediately derive the behaviour of the 
monopole mode |lfj|-[l2]| by computing the evolution of 
the square radius: 



IV. QUADRUPOLE OSCILLATION AND 
RELAXATION TIME APPROXIMATION 

The purpose of this section is to investigate the 
quadrupole mode of a classical gas as well as its cou- 
pling with the monopole oscillation arising in anisotropic 
traps. In this case the solution of the Boltzmann equa- 
tion exhibits damping and one has to deal explicitly with 
the collisional term. In the presence of anistropy, the 
l z = component of the quadrupole is coupled with the 
monopole and one finally finds the following set of cou- 
pled equations : 



d(xi) 
dt 

d(x2 
dt 



2(X3 







- 2<X4> = 



dt 



2(r.v) 



(7) 



In order to obtain a closed set of equations one also needs 
the following equations 



d(r.v) 
dt 



(v 2 ) -u 2 (r 2 } 



and 



(8) 



(9) 



d(X3 
dt 



, 2uj 2 +uj 2 uj\ I- uj 2 , 

(X5> + — =4 -(Xi) + ± (X2) = 



^ HM) + M^ (xi) + !fi+M te) = o 



*») ,M^ te) + M±M fel) = te4oll) 



dt 



(10) 
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where we have defined the quantities 

Xi =r 2 
X 2 = 2z 2 - r\ 
X3 = r ■ v 

X4 = 2zv z - r 



X5 
X6 



2w 2 



(11) 



If the trap is isotropic = u> z ), ( |l0| ) decouples in 
two subsystems. One subsystem refers to the undamped 
monopole oscillations discussed in the previous section. 
The other corresponds to the damped quadrupole mode. 
Notice that collisions affects only the last equation of 
([Tof). Actually only the variable xe = 2v 2 — v\ is not 
a conserved quantity and hence does not satisfy the cri- 
terium (|^) . The collisional contribution can be estimated 
by treating the collisional integral in the so-called relax- 
ation time approximation |Tc[| : 



f-h 



(12) 



where r is a relaxation time. As we will see in the next 
section this time can be explicitly calculated in terms of 
the scattering cross section. The function /^. e . entering 
( |l2| ) is the local equilibrium distribution function and is 
characterized by the same density, velocity field and local 
temperature as the distribution function /. Differently 
from / it has, however, an isotropic velocity distribution. 
As a consequence the average of \6 on /^. e . exactly van- 
ishes and one can write 



coiw 



<X6> 



(13) 



The system ([To]) is now linear, and one can search for 
solutions of the type e tuJt and deduce the associated de- 
terminant. The result is: 



(w 2 -4j 2 )(w 2 -4wi) 



— I cj 4 - -uj' z (buj\ + Acu'i) + 8uj'iuj'i ) = 



(14) 



The first term of (|l4| ) corresponds to the dispersion law 
for the pure collisionlcss regime (lo z t — > oo). In this 
case the eigenfrequencies coincide with the ones predicted 
by the single particle harmonic oscillator hamiltonian: 
locl = 2 u> z and uich = 2 u± . On the contrary, the term 
multiplying 1 /r refers to the pure hydrodynamical regime 
Kt— >0). 

For the spherical trap, one gets whd = v2wq and 2wo 
for the quadrupole and monopole modes respectively. For 
a cigar shaped configuration (A 2 -C 1) the two hydrody- 
namic solutions have the form whd 
•y/10/3 cj± @, while for a disk trap (A 2 
^hd = \/8/3 uj z and \/3 loj_ . 



yl2/5 uj z and 
3> 1) one finds 



Formula ( ^4| ) , which provides the proper interpolation 
between the collisionless and hydrodynamic regimes, can 
be simplified in the case of a spherical, cig ar, and disk 
shaped trap. In fact, the dispersion law ( |14|) can be writ- 
ten in all these limiting cases in the useful form: 



^CL 



7 HD 



CL 



1 + iujf 



(15) 



typical of relaxation phenomena |7|,[l3||. The time f is 
related to r by a simple numerical factor. For example 
f = r for the quadrupole mode in the spherical case, 
and f = 6t/5 for the lowest mode of the cigar shaped 
configuration. A relevant feature of (|l^) is the presence 
of an imaginary part, associated with the damping of 
oscillations. By writing the frequency as lu = uj r + iT one 
finds, assuming r -C u> r , 



r 



2 1 + (u> r f) 2 



(16) 



Notice that the damping depends crucially on the differ- 
ence between the frequencies calculated in the collisional 
and hydrodynamic regimes and exactly vanishes when 
these frequencies coincide. This happens, for example, 
in the monopole case for isotropic trapping, as discussed 
in the previous section |14|. 

In the hydrodynamic limit (u> r f — * 0) the damping 
predicted by (16) takes the form 



Thd 



2^. 



2 

CL 



while in the opposite regime (r 
1 



oo ) one gets 



r 



CL 



2w 



2 

CL 



(17) 



(18) 



CL' 



A maximum for Y is found at uj r f ~ 1, leading to T 



CL 



HDi 



Around this value, the approximation 
leading to (h6|) is no longer accurate, and one should 
rather use (|14|) or jlq) . 



V. GAUSSIAN ANSATZ 

One expects the relaxation time r to be of the same order 
of the inverse of the collision rate 7 co n [fH : 



Tcoll 



n(0)v t h(To 



(19) 



giving the number of collisions undergone by a given 
atom per unit of time. In (|l^) i>th = y/Wo/itm is the 
thermal velocity and ?^(0) the central density. 

An explicit evaluation of this link can be obtained by 
making a gaussian ansatz for the distribution function, 
which allows for an analytic evaluation of the collisional 
contribution to the equations of motion. In fact, a natu- 
ral assumption to describe the l z = modes characterized 
by axial symmetry is given by the parametrization 



3 



/(r,v,i) 



N 
\2tt 



,3 1/2 



with r given by 



e -mUl/2B ±e - m U 2 j28z e - m ( 1± rl+ 7 zZ 2 )/2 



(20) 



where r± = \J x 2 + y 2 and U = v — (v) . Eq. ( f20| ) pro- 
vides a natural generalization of the local equilibrium 
distribution, by introducing a deformation not only in 
coordinate space (taken into account by the 7 param- 
eters), but also in velocity space. The gaussian ansatz 
( pp| ) can be shown to describe exactly the monopole and 
quadrupole oscillations both in the hydrodynamic and 
collisionless regimes. In the former case, the velocity dis- 
tribution is isotropic and hence 9± = Z . For the colli- 
sionless case, an exact solution of the equations of motion 
can be found in the form (|20|). However, in this case, 9± 
is different from 9 Z , corresponding to configurations far 
from local equilibrium. Only in the presence of isotropic 
trapping and for the monopole oscillation the hydrody- 
namic and collisionless solutions coincide. In this case 
the ansatz (^), with 6± — 9 Z and j± — j z provides an 
exact solution of the Boltzmann equation. 

In the limit of small oscillations around the equilibrium 
configuration, the temperature in axial and transverse 
axis can be expanded around the equilibrium tempera- 
ture 9q : 



■X = co 
% = #0 



One then finds: 



(Xe) = -(60 \ 
m 



(21) 



(22) 



showing that (xq) is, as expected, directly sensitive to the 
anisotropy of the velocity distribution. The collisional 
contribution to the equation for xe (last eq. of the system 
(|l0|)), can be also expressed in terms of this anisotropy. 
By inserting eq. (|2^) into eq. (|^) with x = X6 we obtain, 
after linearization, the expression 



(Xe-fcoll) 



/o(2)A Xe 



1 mcr 

32^ m 2 



cPrcPUKpUzlU! - U 2 |d 2 fi/ (1) 



±[m) v + (ui) v -(ui) 1 -(uih)+ 



(u 2 z h - (ul), 



(23) 



where fo is the gaussian ( pOj ) evaluated at equilibium, 
and Axe = (xe)i + (xeh -Xxb)i' - (xeV- After some 
length but straightforward algebra (see Appendix) the 
integral (E3T) can be written in the useful form: 



2 (59, - I 

(Xeicou) = 

t m 



-) (%6> 



(24) 



47coii 



(25) 



Eq. (|2Jj) provides the explicit link between 7 co n and the 
relaxation time for the quadrupole l z — mode. Notice 
that this relationship, provided by the gaussian ansatz 
(pp|), is independent of the deformation of the trap. In 
the next section we will compare this prediction with 
the result of a numerical calculation based on molecular 
dynamics. 



VI. NUMERICAL SIMULATION 

In this section, we present results for the dispersion law 
arising from a numerical analysis, based on a molecular 
dynamics simulation. We consider N — 2 x 10 4 parti- 
cles moving in the potential (j^) . Binary elastic collisions 
are taken into account using a boxing technique |]l6| , |l7f . 
At each time step St, the position of each particle is dis- 
cretized on a square lattice with a step £. The volume 
£ 3 of a box is chosen such that the average occupation 
p occ of any box is much smaller than 1. Collisions occur 
only between two particles occupying the same box, and 
the time step St is adjusted in such a way that the prob- 
ability peon of a collisional event during St is also much 
smaller than 1. We choose typically p occ ~ p co u ~ 5%. 

Initial conditions for exciting the lowest energetical 
mode are obtained by a deformation in coordinate and 
velocity spaces of the cloud along the weak axes, keep- 
ing the phase space density constant. We have checked 
that this method leads to the excitation of only the low- 
est frequency mode. Then, we let the cloud evolve. 
The damped oscillation of the dynamical variable X2 = 
2z 2 — r\ is analysed for different choices of the collision 
rate. As an example, the oscillation frequency and the 
damping for a cigar shaped trap (A = 1/10) are plotted 
(solid circles) respectively on figures (1) and (2). One ob- 
serves that the frequency decreases as the ratio u) z /j co \i 
decreases and tends asymptotically to the hydrodynamic 
value y/12/5 oj z = 1.55 u> z . For large value of uj z Hco\\i 
the frquency instead approaches the collisionless value 
2uj z . By performing a least square fit with formula jl4|), 
we obtain: 



r = (1.31 ±0.07)- 



7coll 



(26) 



which well agree with the gaussian prediction (25). The 
full line on figure (1) and (2) corresponds to eq. (|14|) with 
r given by the gaussian ansatz prediction (|25|). We have 
checked that, for deformed cigar traps (A <C 1), result 
( p6| ) is independent of the value of A. 

We have instead found that, as the deformation param- 
eter A becomes larger, the results of the numerical simu- 
lation exhibits an increasing deviation from the gaussian 
prediction. A discrepancy of the order of 25 % has been 
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observed for A = v°- The physical origin of this devia- 
tion remains to be undestood. In all cases, the gaussian 
ansatz underestimates the coefficient of proportionality 
between the relaxation time r and the inverse of the col- 
lision rate l/7 co n- 



VII. CONCLUSION 

In this work, we have presented an investigation of 
the collective frequencies of a classical gas trapped in 
an harmonic potential well. Starting from the classi- 
cal Boltzmann equation, we have derived a set of cou- 
pled equations for the averages of the relevant dynamic 
variables associated with the monopole and quadrupole 
modes. The collisional term has been treated in the time 
relaxation approximation. The corresponding relaxation 
time was evaluated by a gaussian ansatz for the distri- 
bution function. The quality of the gaussian prediction 
has been checked numerically by a numerical simulation 
based on molecular dynamics. 



APPENDIX A: COLLISIONAL INTEGRAL 

The appendix is devoted to the explicit calculation of 
the collisional integral (|23|). Let us introduce the center 
of mass velocity C and the relative velocity before (V) 
and after (V) collision : 

Ux = C + V/2 
U 2 = C - V/2 
XJ V = C + V/2 

U 2 < = C - V/2 (Al) 

The conservation of kinetic energy during an elastic col- 
lision ensures 



V 2 = V' 2 , 



(A2) 



so that the collisional integral can be rewritten in the 
form 



(Xeicoii) = -(66. - S6±) 



3 TflO Q 

128-kW 2 



<fr d A V d A C d'Sl V /o(l) /o(2) [V z z - V z iY ■ (A3) 



Let us first calculate the angular integral : 

I 2 



Jdn[v?-v?y 



(A4) 



We introduce a reference frame 5ft linked to V, the vec- 
tors of the associated orthonormal basis beeing (a, b, c). 
Without loss of generality we choose a such that V = 
V k, and the z axis in the plane generated by (a, b). The 



relative velocity V is characterized in 5ft by two spherical 
angles (9', ip') : 



V.a=Fcos0' 
V.b = Vsm9' cos <p' 
V.c = V sin 9' sin if' 



(A5) 
(A6) 
(A7) 



Thus 



V' z = V' • z = V z cos 6' + V sin 9' cos ip'(b ■ z) , 

where z is the unit vector of the z axis and 

V' 2 = V 2 cos 2 9' + V 2 sin 2 9' cos 2 (p'(b ■ i) 2 
+ 2VV Z cos 9' sin 9' cos <p'(h ■ z) . 

With our choice of coordinate, one has 

z • z = 1 = (a - z) 2 + (b ■ z) 2 , 

which implies 

^ 2 (b-z) 2 = ^ 2 -i4 2 

and, finally, 

V 2 - V' 2 = V 2 (l - cos 2 9') - sin 2 9' cos 2 V > ' (V 2 - V 2 ) 
2 V V z sin 6' cos 6' cos y' (b • z) . 

By integrating the square of the previous expression, 
one finds 



32-7T /15 A 3 A 5 o r 

In = — V 4 + -V 4 V V 

Q 15 I 8 2 8 4 2 



(A8) 



The calculation of the collisional integral ( |A3| ) is now 
straightforward, and finally yields the result : 



(Xe/coii) = -{56 z - S9 ± )—v th a n(0) , (A9) 
5m 



where Vth = \i 8 flp / 7 rw is the thermal velocity of a par- 
ticle of the gas. (|A9| ) permits to derive ( p4] , p5| ) with 7 co n 
defined by eq. (|f 9|). 
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FIG. 1. Real part of frequency of the l z = mode of a 
classical gas confined in a cigar shaped trap (A = 1/10), ver- 
sus to z /t coii . The solid curve represents the prediction of the 
gaussian ansatz. The circles are the numerical results ob- 
tained with a moleclular dynamics simulation. 
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FIG. 2. Damping of the l z = mode of a classical gas 
confined in a cigar shaped trap (A = 1/10), versus lj z /~{ co \\. 
The notations for the line and the markers are the same as in 
figure 1. 
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